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1.0  SUMMARY 


Episodes  of  hypoxia-like  events  experienced  by  F-22  pilots  after  training  missions  have  raised 
questions  about  pilot  exposures  and  their  effective  mission  accomplishment.  These  episodes 
identified  data  gaps  in  the  understanding  of  pilot  exposures  in  high  performance  aircraft.  We, 
therefore,  need  to  understand  how  diminished  levels  of  oxygen  at  lower  than  atmospheric 
pressure  may  cause  physiological  hypoxia,  and  how  temporary  loss  of  this  metabolic  and 
physiologic  requirement  relates  to  alterations  in  cognition. 

Responses  to  hypoxia  take  place  at  multiple  levels  within  the  human  body.  The  initial  response  is 
through  the  carotid  body  altering  physiology  (breathing  rate,  regional  blood  flow).  Tissue  level 
responses  involve  regulating  demand  for  oxygen  through  altered  function.  Key  to  understanding 
these  responses  both  at  the  physiological  and  tissue  (brain)  level  are  oxygen  sensing  pathways, 
particularly  those  involving  the  oxygen  sensing  cytosolic  gene  transcription  factor  HIF-la.  These 
linked  oxygen  sensing  pathways  need  to  be  understood  mechanistically  and  quantitatively  for 
defining  response  changes  and  assessing  changes  in  cognition. 

A  mathematical  model  was  developed  to  integrate  these  experimental  measurements  and  provide 
a  description  of  the  mechanisms  of  hypoxia  and  the  mammalian  organism’s  response.  Our  model 
focuses  on  oxygen  delivery  to  the  brain.  It  includes  delivery  of  oxygen  to  the  blood  via 
hemoglobin  binding,  breathing  rate  and  cardiac  output,  delivery  of  oxygen  from  the  blood  to  the 
brain  tissue,  and  finally  the  response  of  the  HIF-la  signaling  pathway  in  terms  of  alterations  in 
gene  expression  and  other  end-points.  The  model  simulates  time  courses  of  tissue  (brain) 
oxygen,  intracellular  HIF-la,  transport  to  the  nucleus,  binding  to  hypoxia  responsive  element 
(HRE)  and  generation  of  transcription  signal,  during  altered  blood  oxygen  levels.  Model 
development  is  an  iterative  process;  the  model  described  here  is  sufficient  to  describe  these 
underlying  mechanisms,  and  provides  a  firm  basis  for  further  development,  should  the  need  arise 
to  incorporate  additional  data.  Several  such  potential  data  extensions  that  would  advance  the 
model  are  outlined. 

The  HIF-la  model  is  partially  parameterized  based  on  experimental  data  in  the  literature,  as  well 
as  existing  models  related  to  or  including  the  HIF-la  signaling  pathway.  The  model  has  been 
fitted  to  (literature)  steady- state  dose/response  data  for  vascular  endothelial  growth  factor 
(VEGF)  and  VEGF  mRNA.  The  model  also  describes  angiogenesis  data  (longer  term  changes  in 
capillary  densities  in  response  to  prolonged  hypoxia  exposures),  while  assuming  linearity  of 
dose/response  for  gene  transcription  process. 

In  vitro  cultured  cell  studies  are  useful  for  elucidating  elements  of  these  processes.  These  data 
need  to  be  integrated  and  applied  (extrapolated)  to  relevant  human  exposure  scenarios  in  order  to 
assess  the  relevance  of  hypoxia  to  the  abovementioned  Air  Force  (AF)  issues.  Mathematical 
models  of  the  kind  described  here  are  essential  both  for  the  understanding  of  in  vitro  results 
investigating  underlying  mechanisms,  and  in  order  to  extrapolate  these  results  to  in  vivo  and 
relevant  human  scenarios. 
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This  report  continues  and  concludes  the  work  initiated  in  a  previous  Technical  Report  AFRL- 
RH-WP-TR-  2016-0086:  Host  -  HIF-la  Pathway  and  Hypoxia:  In  Vitro  Studies  and 
Mathematical  Model ,  by  Robinson  and  others. 
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2.0  INTRODUCTION 


2.1  Background 

There  are  currently  issues  with  multiple  aircraft  in  the  US  armed  forces,  including  the  AF  F-22, 
F-35,  A-10  and  Navy  F/A-18,  possibly  related  to  hypoxia,  as  well  as  other  conditions,  such  as 
hyperoxia  and  hypobaria,  that  may  negatively  impact  pilot  performance  (Martin  et  al.,  2012).  In 
addition,  pilots  of  the  high-altitude  U-2  reconnaissance  aircraft  may  suffer  brain  lesions  (white 
matter  hyperintensities)  and  cognitive  deficits  (McGuire  et  al.,  2014).  In  light  of  these  issues,  we 
need  to  understand  how  oxygen  delivery  (as  a  result  of  both  concentration  and  air  pressure 
changes)  relates  to  alterations  in  brain  functioning.  In  addition  to  oxygen  alone,  we  also  need  to 
understand  the  kinetics  of  potential  contaminants  under  different  oxygen  conditions. 
Physiologically-based  pharmacokinetic  (PBPK)  models  have  been  developed  in  our  lab  for 
numerous  chemicals  and  mixtures  of  concern.  However,  no  PBPK  models  of  combined  oxygen 
and  chemical  contaminant  exposure  has  been  developed. 

Organisms  respond  to  low  oxygen  levels  in  the  air  in  a  number  of  ways.  An  initial  response  is 
through  carotid  body  (CB)  that  mediates  alterations  of  physiology  (breathing,  heart  rate/blood 
flow).  Tissue  level  responses  involve  regulating  demand  for  oxygen  through  altered  function. 

Key  to  understanding  these  responses  both  at  the  physiological  and  tissue  (brain)  level  are 
oxygen  sensing  and  response  pathways,  particularly  those  involving  hypoxia  inducible  factor- 1- 
alpha  (HIF-la).  These  linked  pathways  need  to  be  understood  mechanistically  and  quantitatively. 
In  vitro  studies  are  useful  for  elucidating  elements  of  these  processes.  These  approaches  need  to 
be  integrated  and  applied  (extrapolated)  to  relevant  human  exposure  scenarios  in  order  to  assess 
the  relevance  of  hypoxia  to  the  abovementioned  AF  issues.  Initial  in  vitro  studies  have  been 
described  in  a  previous  Technical  Report  AFRL-RH-WP-TR-2016-0086:  Host  -  HIF-la 
Pathway  and  Hypoxia:  In  Vitro  Studies  and  Mathematical  Model,  by  Robinson  et  al. 

Compensatory  mechanisms  for  hypoxia  in  brain  (and  other  tissues)  work  at  a  number  of  time 
scales,  and  include  physiological  and  cellular  responses.  Physiological  responses  are  fast- acting, 
and  prevent  cellular  hypoxia  (within  limits).  These  responses  are  typically  mediated  by  the 
carotid  (and  aortic)  body,  which  is  a  sensory  organ  for  detecting  arterial  blood  O2  levels  and 
mediating  systemic  cardiac,  vascular  and  respiratory  responses  to  hypoxia.  Carotid  body  (CB) 
response  travels  via  afferent  nerves  that  carry  signals  back  from  the  carotid  and  aortic  bodies  to 
the  brainstem,  which  responds  accordingly  by,  for  example,  increasing  the  ventilation  rate. 

Figure  1  illustrates  this  adaptive  response  to  hypoxia. 

Responsiveness  of  the  CB  to  acute  hypoxia  relies  on  the  inhibition  of  02-sensitive  K+  channels  in 
glomus  cells,  which  leads  to  cell  depolarization,  Ca2+  entry  and  release  of  transmitters  that 
activate  afferent  nerve  fibers;  the  molecular  mechanisms  underlying  K+  channel  modulation  by 
O2  tension  may  involve  disruption  of  the  prolylhydroxylase/  hypoxia  inducible  factor  pathway 
(Lopez-Bameo  et  al.,  2008).  Evidence  that  HIF-la  is  involved  with  carotid  body  responses  to 
hypoxia  includes  studies  with  HIF-la  +/-  mice  (Kline  et  al.  2002;  Peng  et  al.  2006),  in  which  the 
hypoxia  response  was  severely  impaired.  However,  there  is  some  disagreement,  with  some 
researchers  thinking  it  unlikely  that  acute  responses  to  hypoxia  by  the  carotid  body  (such  as 
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regulation  of  breathing  and  heart  rate)  are  mediated  directly  by  the  HIF-la  pathway  (Wyatt, 
personal  comm.).  In  the  case  of  chronic  hyperoxia,  the  02-sensitive  response  is  diminished 
(leading  to  cell  hyperpolarization)  while  the  hypercapnic  (blood  CO2  retention)  response  is 
augmented  (Sinski  et  al.,  2014). 


HIF1,2a  models, 
applied  to  suitable 
in  vitro  systems 


Figure  1.  Schematic  of  the  Physiological  Response  to  Hypoxia  via  the  Carotid  Body.  Figure 
also  shows  the  tissue  response  to  hypoxia,  both  mediated  by  the  HIF-la  signaling  pathway 


When  such  first  tier,  physiological  protection  mechanisms  are  inadequate,  cellular  hypoxia 
ensues,  and  affected  cells  such  as  neurons  adapt  to  the  hypoxic  conditions  by  altering 
homeostasis  through  gene-regulated  cellular  processes  such  as  metabolism,  apoptosis, 
angiogenesis  and  cell  proliferation  and  differentiation. 

We  propose  to  develop  a  mechanism-based  model  for  an  organism’s  response  to  hypoxia.  Our 
model  involves  two  interlinked  parts:  Firstly,  a  whole  body  oxygen  delivery  model  that 
determines  blood  oxygen  levels  from  inhaled  air;  and  secondly,  a  HIF-la  signaling  model  that 
determines  the  tissue  response. 
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2.2  Oxygen  Transport  to  Tissues 


About  98.5  percent  of  the  oxygen  in  the  arterial  blood  in  a  healthy  human  at  sea-level  pressure  is 
bound  to  hemoglobin,  while  the  remaining  1.5  percent  is  dissolved  in  the  plasma  and  other  blood 
liquids.  The  hemoglobin  molecule  is  the  primary  transporter  of  oxygen  in  mammals  and  many 
other  species.  Hemoglobin  has  an  oxygen  binding  capacity  of  1.368  ±  0.017  (mean  ±  standard 
deviation)  mL  O2  per  gram  hemoglobin  (Dijkhuizen  et  al.,  1977).  Table  1  shows  Normal  values 
of  oxygen  concentration  in  healthy  humans,  expressed  in  mmHg  (partial  pressure),  percent  O2, 
and  mmol/L. 


Table  1:  Normal  (Mean)  Values  of  Oxygen  Concentration  in  Healthy  Humans 


Location 

p02  (mmHg) 

%02 

mmol/L 

Air 

160 

21.0 

8.3 

Trachea 

150 

19.7 

7.8 

Alveoli:  100-120  mmHg 

110 

14.4 

5.7 

Lung  inflow  (arterial):  80-100mmHg 

90 

11.8 

4.7 

Lung  outflow  (venous) 

40 

5.2 

2.1 

Arterial 

100 

13.1 

5.2 

Venous 

40 

5.2 

2.1 

Brain# 

35 

4.6 

1.8 

Notes:  Adapted  from  Carreau  et  al.  (201 1);  #  Dings  et  al.  (1998)  measured  brain  oxygen 
concentrations  as  a  function  of  depth  using  a  brain  inserted  probe  in  human  patients,  and 
observed  mean  p02  of  23.8  ±8.1  mmHg  at  22  to  27  mm  below  the  dura,  25.7  ±  8.3  mmHg  at  17 
to  22  mm,  33.0  ±  13.3  mmHg  at  12  to  17  mm  and  33.3  ±  13.3  mmHg  at  7  to  12  mm.  The 
“scientific  community  agrees  to  say  that  if  the  brain  tissue  p02  is  above  35  mmHg  (4.6  percent), 
normal  oxygenation  of  the  brain  tissue  should  be  assured”  (Carreau  et  al.,  201 1). 


As  an  example  of  brain  tissue  and  blood  oxygen  rapid  responses  to  short  periods  of  hypoxia, 
Figure  2  shows  pCfi  levels  in  arterioles,  venules,  and  tissue  sites  in  the  venous  capillary  network 
during  a  1 -minute  control  period,  1 -minute  reduction  of  inspired  oxygen  fraction  from  21  percent 
to  7  percent,  and  2-minute  recovery  period  (Johnson  et  al.,  2005). 
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Figure  2.  pCh  Levels  in  Arterioles,  Venules,  and  Tissue  Sites  in  the  Venous  Capillary 
Network.  Figure  shows  response  during  a  1 -minute  control  period,  1 -minute  reduction  of 
inspired  oxygen  fraction  from  21  percent  to  7  percent,  and  2-minute  recovery  period  at  0.21. 
Data  points  are  means  +/-  SE.  Data  from  Johnson  et  al.  (2005). 


The  literature  is  rich  with  models  of  oxygen  transport  and  cardiopulmonary  control,  at  various 
levels  of  physiological  detail  (Grodins  et  al.,  1967;  Duffin  et  al.,  2000;  Topor  et  al.,  2004;  Ward, 
2007;  Wolf  and  Garner,  2007;  Benson  et  al.,  2013;  Cheng  et  al.,  2016).  For  example,  Topor  et 
al.  (2004)  have  developed  a  mathematical  model  that  describes  multiple  mechanisms  associated 
with  control  of  breathing  in  response  to  hypoxia.  Although  these  models  are  adequate  for 
predicting  oxygen  transport,  the  model  equations  are  in  quite  a  different  format  than  those  of  our 
typical  PBPK  models  for  chemical  contaminants,  making  the  combination  of  the  models 
difficult.  Therefore,  a  goal  of  this  effort  was  to  investigate  the  applicability  of  our  PBPK 
approach  to  the  transport  of  oxygen. 

The  impact  of  chemical  contaminants  on  oxygen  transport  is  often  through  an  allosteric  shift  in 
oxygen-hemoglobin  binding.  For  example,  organic  phosphates,  such  as  2,3-Diphosphoglycerate 
(2,3-DPG)  binds  to  hemoglobin  and  rearranges  it  into  a  tense  state,  which  decreases  its  affinity 
for  oxygen.  Carbon  monoxide  (CO)  is  also  a  classic  example  of  a  compound  interfering  with  the 
oxygen  transport  function  of  the  blood  by  combining  with  hemoglobin  to  form  hemoglobin 
(COHb)  (see,  for  example,  Berne  et  al.,  2004).  CO  has  approximately  200  times  the  affinity  for 
hemoglobin  that  oxygen  does  and  for  that  reason  even  small  amounts  of  CO  can  tie  up  a  large 
proportion  of  the  hemoglobin  in  the  blood  making  it  unavailable  for  oxygen  carriage.  If  this 
happens,  the  p02  of  the  blood  and  the  hemoglobin  concentration  will  be  at  normal  levels,  but  the 
the  oxygen  concentration  and  transport  to  the  tissues  will  be  grossly  reduced. 
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2.3  HIF-la  Pathway 


In  both  the  carotid  bodies  (Prabhakar,  2013),  and  other  tissues,  the  primary  sensing  molecule  is 
the  hypoxia  inducible  factor  (HIF-la)  transcription  factor  protein.  The  presence  of  oxygen 
activates  prolyl  hydroxylases  (PHD)  to  hydroxylate  HIF-la,  leading  to  HIF-la  protein 
degradation.  Under  low  oxygen  conditions,  hydroxylase  inactivity  allows  HIF-la  protein  to 
accumulate,  which  in  turn  initiates  target  gene  expression.  HIF-la  dimerizes  with  HIF-ip  and 
translocates  to  the  nucleus  where  it  associates  with  CBP  (CREB  binding  protein)  and  DNA 
polymerase  II  and  then  binds  DNA  to  initiate  transcription,  up-regulating  expression  of  apoptotic 
cell  death  pathway  genes,  angiogenic  genes,  erythropoeitic  genes  and  other  protein  genes. 

Figure  3  shows  the  interactions  of  the  HIF-la  protein  under  conditions  of  normoxia  and  hypoxia. 
Under  normoxic  conditions,  hydroxylation  of  Pro-402  and  Pro-564  by  the  HIF-la  prolyl 
hydroxylases  (PHD)  allows  binding  of  von  Hippel-Lindau  tumor  suppressor  protein  (VHL), 
which  permits  HIF-la  ubiquitination  and  degradation,  maintaining  low  levels  of  HIF-la  in  the 
cytosol.  At  the  same  time,  hydroxylation  of  Asn-803  by  FIH1  (factor  inhibiting  HIF-1)  prevents 
the  interaction  of  HIF-la  with  the  cofactors  p300  and  CBP  (CREB  binding  protein).  Under 
hypoxic  conditions,  hydroxylation  is  less  likely,  VHL  does  not  bind  and  HIF-la  builds  up  in  the 
cytosol.  At  the  same  time,  binding  of  p300/CBP  cofactors  facilitates  HIF-la  protein 
translocation  to  the  nucleus  and  subsequent  gene  transcription  (Semenza,  2004). 

A  number  of  experimental  studies  can  be  done  to  elucidate  the  mechanisms  associated  with  both 
the  physiological  and  cellular  responses  to  hypoxia.  In  vivo  studies  can  be  used  to  explore  the 
physiological  responses  mediated  by  the  carotid  body.  In  addition,  there  is  the  possibility  of 
looking  at  HIF  signaling  in  cultured  Type  1  carotid  bodies.  In  general,  in  vitro  studies  can  be 
used  to  study  hypoxia  responses  in  specific  cell  types  such  as  neurons.  For  example,  high- 
throughput  in  vitro  screening  can  be  used  to  monitor  expression  of  HIF-la,  determine  proteomic 
profiles  of  hypoxic  events  and  identify  potential  protein  markers  for  further  investigation  from 
brain  tissues  collected  from  in  vivo  studies.  Immunohistochemical  staining  (IHC)  of  whole  brain 
from  altitude  studies  was  performed  to  identify  regions  most  affected  by  hypoxia.  A  HIF-la 
specific  antibody  can  be  used  to  see  cytosolic  and  nuclear  accumulation  of  HIF-la.  Protein 
biomarkers  identified  from  in  vitro  studies  can  also  be  evaluated  by  IHC.  Finally,  it  is  also 
possible  to  measure  mRNA  levels  of  relevant  proteins  in  the  hypoxia  signaling  pathways. 
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Normoxia: 


Pro-402  | 

I  Pro-564  | 

1  Asn-803  1 

HIFIa  Protein 

TAD  N 

TAD  C 

1 

HIFIa  Degradation 


Hypoxia: 


i 

Gene  Transcription 


1 

Angiogenesis  etc. 


Figure  3.  Binding  Sites  on  HIF-la  Protein.  N-  and  C-terminal  transactivation  domains  (TAD 
N  and  TAD  C)  are  regions  for  HIF-la  degradation  and  transcriptional  activation,  respectively. 
Under  normoxic  conditions,  hydroxylation  of  Pro-402  and  Pro-564  by  the  HIF-la  prolyl 
hydroxylases  (PHD)  allows  the  von  Hippel-Lindau  tumor  suppressor  protein  (VHL)  to  bind, 
which  permits  HIF-la  ubiquitination  and  degradation,  maintaining  low  levels  of  HIF-la  in  the 
cytosol.  At  the  same  time,  hydroxylation  of  Asn-803  by  FIH1  (factor  inhibiting  HIF-1)  prevents 
the  interaction  of  HIF-  la  with  the  cofactors  p300  and  CBP  (CREB  binding  protein).  Under 
hypoxic  conditions,  hydroxylation  is  less  likely,  VHL  does  not  bind  and  HIF-la  builds  up  in  the 
cytosol.  Also,  binding  of  p300/CBP  cofactors  facilitates  translocation  to  the  nucleus  and 
subsequent  gene  transcription.  (Adapted  from  Semenza,  2004). 
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2.4  Summary  of  Existing  Mathematical  Models 


Table  2  shows  a  summary  of  some  existing  models  for  the  HIF-la  pathway.  At  this  point,  our 
current  model  does  not  include  hydroxylation  of  Asn-803  by  FIH1  (factor  inhibiting  HIF-1), 
which  prevents  the  interaction  of  HIF-la  with  the  cofactors  p300  and  CBP  (CREB  binding 
protein).  Under  hypoxic  conditions,  hydroxylation  is  less  likely,  VHL  does  not  bind  and  HIF-la 
builds  up  in  the  cytosol.  At  the  same  time,  binding  of  p300/CBP  cofactors  facilitates 
translocation  to  the  nucleus  and  subsequent  gene  transcription  (Semenza,  2004)  -  see  Figure  3.  A 
number  of  models  do  include  this  process  (see  Table  2).  However,  the  model  described  in  the 
current  Report  is  sufficient  to  describe  the  build-up  of  cytosolic  HIF-la  under  hypoxic 
conditions,  and  provides  a  firm  basis  for  further  model  development,  should  the  need  arise  to 
incorporate  or  describe  additional  data  (or  as  such  data  become  available). 


Table  2.  Brief  Summary  of  Some  Existing  Models  of  the  HIF-la  Pathway 


Model  Reference 

Model  Description 

Notes 

Kohn  et  al.,  2004 

HRE  occupancy  and/or  mRNA  expression  in  response 

Initiated  mathematical  modelling  of  the 

to  changes  in  oxygen  concentration. 

HIF  pathway. 

Kooner  et  at,  2006 

HRE  occupancy  and/or  mRNA  expression  in  response 
to  changes  in  oxygen  concentration. 

Extended  Kohn's  model  with  role  of 
oxygen  in  mediating  reactions,  as  well  as 
its  transfer  to  nucleus. 

Yu  et  at,  2007 

HRE  occupancy  and/or  mRNA  expression  in  response 

Uses  "extreme  pathway  analysis"  rather 

to  changes  in  oxygen  concentration. 

than  solving  system  of  ODEs 

Qutub  &  Popel  2006 

Sensitivity  of  HIF  protein  to  hydroxylation  cofactors 
and  PHD 

Qutub  &  Popel  2007 

Temporal  effect  on  HIF  protein  stabilisation  by 
succinate  inhibition  and  PHD  negative  feedback. 

Yucel  &  Kumaz  2007 

Sensitivity  of  the  angiogenic  behaviour  of  a  cancer  cell 
to  PHD  and  FIH. 

Dayan  et  at,  2009 

FIH  controls  a  switch  between  C-TAD  and  N-TAD  HIF 
target  gene  repertoires 

Schmierer  et  at,  2010 

Sequestration  of  FIH  modulates  FIH  activity 

Nguyen  et  at,  2013 

Global  temporal  dynamics  of  the 

Estimation  of  some  parameters  to  in- 

HIF:  PHD:  FIH  network. 

house  experimental  data 

Adapted  from  Cavadas  et  ah,  2013 


3.0  METHODS 


3.1  Oxygen  Delivery  Model 

Table  1  above  shows  the  decrease  in  oxygen  concentration  in  the  air  as  it  passes  from  the 
external  environment,  down  the  trachea  into  the  alveolar  region.  This  decrease  is  primarily  due  to 
the  increase  in  water  vapor  concentration  in  the  air  as  it  becomes  progressively  more  humidified. 
Oxygen  in  the  lung  alveolus  then  partitions  into  the  blood  and  binds  to  hemoglobin.  Figure  4 
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shows  the  saturation  curve  for  hemoglobin  (A),  together  with  the  (computed)  relationship 
between  free  and  bound  oxygen  at  equilibrium  (B).  In  (A),  the  solid  line  shows  the  degree  of 
hemoglobin  saturation  (left  y-axis)  and  the  concentration  of  the  bound  component  in  the  blood 
(right  y-axis),  as  a  function  of  pCh.  The  broken  lines  (with  the  right-hand  axis)  show  the 
contribution  of  the  unbound  (free,  dissolved)  component  both  as  it  contributes  to  the  total  (upper 
curve)  and  its  absolute  concentration  (lower  curve). 


(A) 


Total  02 


Po2  mmHg 


(B) 


Free  Oxygen  (ml_/L) 


Figure  4.  Hemoglobin  and  Oxygen  Binding  Curves.  (A)  The  solid  line  shows  the  degree  of 
hemoglobin  saturation  (left  axis)  and  the  concentration  of  the  bound  component  in  the  blood 
(right  axis)  as  a  function  of  p02.  The  broken  lines  (with  the  RH  axis)  show  the  contribution  of 
the  unbound  (free,  dissolved)  component  both  as  it  contributes  to  the  total  (upper  curve)  and  its 
absolute  concentration  (lower  curve)  (adapted  from  Collins  et  al.,  2015).  (B)  Relationship 
between  bound  (y-axis)  and  free  (x-axis)  oxygen  concentrations. 


Figure  5  shows  a  schematic  of  the  oxygen  transport  portion  of  our  model.  Oxygen  enters  the 
alveolus  in  air  saturated  with  water  vapor  at  a  partial  pressure  of  47  mmHg,  at  a  ventilation  rate 
of  Qp.  Oxygen  diffuses  into  the  blood  with  a  relative  affinity  or  partition  coefficient  of  PB.  We 
assume  that  oxygen  in  the  blood  rapidly  binds  to  hemoglobin,  so  that  free  (Cr)  and  bound  (Crs) 
oxygen  in  the  blood  are  in  instantaneous  equilibrium.  The  oxygenated  arterial  blood  (Ca  =  Cf  + 
Cb)  is  carried  from  the  lungs  to  the  rest  of  the  body  by  the  cardiac  output  blood  flow  ( Qc ).  Once 
the  oxygenated  blood  reaches  the  brain,  the  higher  pH  in  that  tissue  (due  to  local  carbon  dioxide 
accumulation  from  metabolism)  facilitates  the  release  of  bound  oxygen  for  diffusion  into  the 
tissue.  In  our  model  we  assume  that  both  bound  and  free  oxygen  are  equally  available  for  uptake, 
which  is  limited  by  an  effective  permeability-surface  area  product  (PA)  and  a  brain  partition 
coefficient  (PBr).  Oxygen’s  dissociation  from  hemoglobin  is  heavily  influenced  by  the  tissue  pH, 
therefore  the  brain’s  pH  is  incorporated  into  PBr.  Inside  the  tissues,  oxygen  is  removed  via  a 
saturable  metabolic  rate  ( VMax ,  KM).  The  blood  from  the  brain  and  rest  of  the  body  is  combined 
in  a  mixed  venous  blood  compartment  (concentration  Cv )  and  returned  to  the  alveolar  region. 

The  model  was  implemented  in  R  (R  Core  Team,  2016),  and  the  code  is  shown  in  Appendix  A. 
Supporting  R  scripts  are  shown  in  Appendix  B. 
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Figure  5.  Schematic  of  the  Oxygen  Transport  Model.  Oxygen  in  the  lung  alveolus  diffuses 
into  the  blood,  binds  to  hemoglobin,  and  is  transported  to  the  tissues  where  it  is  metabolized 
(converted  to  CO2). 


Oxygen  binding  to  hemoglobin  is  described  using  a  Hill  equation  relating  free  oxygen  dissolved 
in  plasma  to  bound  oxygen: 


Bmax*Cpn 

KDn+CFn 


(Equation  1) 


The  parameters,  Bmax,  KD,  and  n,  are  determined  with  the  help  of  the  equation  for  the  total 
oxygen  content  of  the  blood: 


C  —  1.39  *  Hgb  *  Sa02  +  0.03  *  Pa02  (Equation  2) 

The  first  term  gives  the  amount  of  oxygen  bound  to  hemoglobin  (Cs),  where  Hgb  is  the  amount 
of  hemoglobin  ( g/L ),  SaOi  is  the  fraction  saturation  of  Hgb,  and  1.39  is  the  amount  of  oxygen 
per  gram  of  Hgb.  The  second  term  is  the  amount  dissolved  (Cf),  where  0.03  is  the  solubility  of 
oxygen  in  the  blood,  and  Pa02  is  the  partial  pressure  of  oxygen  in  the  blood.  [Note  that  SaOi  is  a 
function  of  PaOi,  which  is  in  fact  the  oxygen  dissociation  curve.] 
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A  series  of  PaOi  were  generated  and  used  in  turn  to  generate  a  series  of  corresponding  Cb  and  Cf 
values,  using  typical  values  for  a  healthy  adult  for  Hgb ,  Sa02,  and  PaOi.  The  values  were  then 
used  to  obtain  parameters  fits  for  the  above  Hill  equation  using  the  “nls()”  function,  a  nonlinear 
regression  function  in  the  software  R. 

Now  what  remains  is  to  describe  how  Cf  was  determined.  The  rate  equation  for  the  arterial  blood 
is  the  following: 

rlA 

-=Qp*(Cj-  CF/PB)  +  Qc  *  (Cv  -  CA)  (Equation  3) 

where  Qp  is  the  pulmonary  blood  flow  rate,  Ci  is  the  inspiratory  alveolar  oxygen  concentration, 
Cf/PB  ( Cx)  is  the  exhaled  oxygen  concentration,  Qc  is  the  cardiac  blood  output,  Cv  is  the  mixed 
venous  blood  concentration,  and  Ca=  Cf  +  Cb  is  the  total  arterial  blood  concentration.  We 
assume  that  equilibration  of  oxygen  between  the  alveoli  and  blood  is  rapid,  and  thus  set  the 
derivative  equal  to  zero,  arriving  at  an  equation  that  is  solved  implicitly  for  Cf  at  each  step  of  the 
integration.  With  Cf  computed,  Cb  can  be  calculated,  and  finally  the  Ca  leaving  the  gas  exchange 
region. 


3.2  HIF-la  Mathematical  Model  for  Brain 

A  mathematical  model  was  developed  to  integrate  these  experimental  measurements  and  provide 
a  description  of  the  mechanisms  of  hypoxia  and  the  mammalian  (or  human)  organism’s  response. 
The  initial  iteration  of  our  model  focuses  on  the  brain.  It  includes  delivery  of  oxygen  to  the  cell, 
and  the  response  of  the  HIF-la  signaling  pathway  in  terms  of  alterations  in  gene  expression  (see 
Figure  6). 

Our  model  is  necessarily  simplified.  Our  aim  is  to  capture  the  essential  behavior  of  the  system  in 
a  way  that  relates  to  the  available  and  relevant  experimental  data.  Our  initial  model  is  partially 
parameterized,  and  initial  estimates  are  made  from  parameter  values  in  order  to  elicit  appropriate 
model  behavior.  As  experimental  data  is  incorporated  into  the  model,  these  parameter  values  will 
be  updated  and  refined  in  an  iterative  manner.  Briefly,  the  model  assumes  that  oxygen  movement 
from  the  blood  into  the  intracellular  space  is  determined  by  a  permeability-surface  area  (PA) 
product.  Once  in  the  cell,  it  is  removed  via  respiration  at  a  maximum  (saturable)  rate  Vresp.  PHD 
hydroxylates  HIF-la,  with  intracellular  levels  of  O2  as  a  co-substrate  (Semenza  et  al.,  2004). 

HIF-la  levels  are  determined  from  a  balance  of  a  constant  zero  order  synthesis  rate  Khif,  and  a 
PHD-dependent  degradation  rate  Kphd.  In  addition,  HIF-la  is  transported  into  the  nucleus,  where 
it  impacts  gene  expression.  We  assume  a  maximum  transport  rate  Tmax,  and  a  half- saturation 
concentration  of  Kt. 
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Figure  6.  Schematic  of  the  HIF-la  Pathway  in  Cells  Such  as  Neurons.  The  pathway  responds 
to  cellular  hypoxia  by  altering  gene  expression. 


The  equations  describing  O2  and  HIF-la  levels  in  the  intracellular  compartment  (volume  V<)  and 
the  nucleus  (volume  14)  are  thus: 


Y  d[o2] c 
c  dt 


P/t.([02]„ 


[02U 


Vresp-[02]c  _  R  \pHD]  Vc.[mF].[Q2]c 

Kresp  +  I02]c  PHD'  J'  ([H/F]+FmH/F)(Fm02  +  [02]c) 


d[HiF]_  rpuni _ [h;f].[o2]c _ rmax[tf/F] 

dt  HIF  PHD J'([H/F]+KmH/F)(Km02+[02]C)  [H/F]+FT 


(Equation  4) 
(Equation  5) 


d[w/F]m 

dt 


TmaxlHIF] 

[HIF]+Kt 


-  K£eg[HIF]n  -  Kbm[HIF]n[HRE]  , 


(Equation  6) 


where  [ 02]P  and  [ 02]c  are  the  concentrations  (partial  pressures)  of  oxygen  in  the  plasma  and  cell 
intracellular  compartments,  [PHD]  and  [HIF]  are  the  intracellular  concentrations  of  prolyl 
hydroxylase  and  HIF-la,  respectively,  and  [HIF\n  is  the  concentration  of  HIF-ndegla  in  the 
nucleus.  Figure  7  shows  a  schematic  of  the  model,  with  the  modeled  components  in  blue,  and 
parameters  or  fixed  values  in  red. 
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Once  in  the  nucleus,  HIF-la  is  assumed  to  bind  with  hypoxia  responsive  element  (HRE)  on  the 
DNA  with  a  bimolecular  rate  constant  Kbm,  and  is  degraded  according  to  a  first  order  rate 
constant  Kndeg  (equation  4).  The  binding  of  HIF-la  to  HRE  initiates  a  transcription  signal  TC, 
which  is  assumed  to  mediate  gene  transcription.  TC  is  assumed  to  spontaneously  decay  at  a  rate 
Ktc: 


—  =  Kbm  [HIF]n  [HRE]  -  Ktc.  TC.  (Equation  7) 

These  equations  were  incorporated  into  Berkeley  Madonna  (see  Robinson  et  al.,  2016)  and  also 
translated  into  “R”,  and  combined  with  the  oxygen  transport  model  (Section  3.1  above)  for  the 
simulations  that  follow  -  see  Appendix  A  and  B. 


Figure  7.  Schematic  Representation  of  the  Kinetic  Model.  Components  whose  concentration 
time  course  is  modeled  in  response  to  changes  in  blood  oxygen  levels  are  shown  in  blue,  while 
parameters  or  fixed  values,  that  ultimately  need  to  be  experimentally  determined,  or  fitted,  are  in 
red. 


3.3  Parameter  Values 

A  number  of  parameters  for  the  model  were  gleaned  from  the  published  literature.  KmHiF,  the 
Michaelis-Menten  constant  for  HIF-la  as  a  substrate  of  PHD  is  reported  as  100  nM  (Hirsila  et 
al.,  2003;  Koivunen  et  al.,  2007;  Kaelin  and  Ratcliffe,  2008,  as  reported  by  Nguyen  et  al.,  2013). 
Kphd,  the  catalytic  rate  constant  for  PHD-mediated  hydroxylation  of  HIF-  la  was  taken  to  be 
0.045  nM-1  s-1  (from  Nguyen  et  al.,  2013).  Km02,  the  Michaelis-Menten  constant  for  O2  as  a 
substrate  of  PHD  is  250  pM  (Hirsila  et  al.,  2003;  Koivunen  et  al.,  2007;  Kaelin  and  Ratcliffe, 
2008,  reported  by  Nguyen  et  al.,  2013).  A  typical  neuron  has  a  volume  of  105  urn2  (diameters 
range  from  4-100  um  in  diameter  -  Davies,  2002),  while  the  volume  of  nucleus  is  assumed  to 
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be  5-25  percent  of  volume  of  neuron  (we  assume  an  average  of  about  15  percent).  The  basal 
HIF-la  synthesis  rate  used  was  0.005  nM  s'1  (Nguyen  et  al.,  2013). 


3.4  Integrating  the  Oxygen  Transport  and  HIF-1  a  Models 

The  oxygen  transport  and  HIF-1  01  models  were  combined  by  embedding  the  HIF-1  01  model  in 
the  brain  tissue  compartment  of  the  oxygen  delivery  model.  (A  number  of  parameters  appear  in 
both  models,  sometimes  with  different  names;  see  Appendix  C:  Model  Parameters  for 
disambiguation).  The  rate  equations  for  the  brain  tissue  and  cytoplasmic  HIF-1  01  were  modified 
as  follows  to  integrate  the  two.  The  equation  for  cellular  oxygen  in  the  HIF-la  model  (equation 
4)  was  modified  to  be  consistent  with  the  parameters  in  the  oxygen  delivery  model,  to  yield: 


dABr,T  _  pA  (py  _  Q  /p  \  _  VMaxBr-CBr,T  _  p  FPHp[PHD][HIF]-CBr,T 

dt  ^  Br  Br,T  Br '  KM+CBr,T  (iHlF]+KmHIF)(CBrT+Krno2) 


(Equation  8) 

where  A  is  amount  (mg)  of  oxygen,  C  is  concentration  (mg/L),  and  the  subscript  Br  indicates 
brain,  and  T  is  for  the  brain  tissue  subcompartment.  F  is  a  conversion  factor  from  nM/second  to 
mg/hour,  making  use  of  the  volume  of  the  brain  tissue.  The  HIF-1  01  rate  equation  was  modified 
by  replacing  [02]c  with  Cb,-.t.  The  remaining  HIF-1  Ot  model  equations  (6  and  7)  were  added 
unchanged  to  the  transport  model. 


4.0  RESULTS 


4.1  Time-Course  Simulations 

The  simulation  of  oxygen  transport,  HIF-la  kinetics,  and  transcription  factor  production  using 
the  combined  model  is  shown  in  Figures  8  and  9.  The  scenario  simulated  is  normoxia,  or  21 
percent  oxygen  inhaled  for  6  hours,  hypoxia,  or  7  percent  oxygen  inhaled  for  4  hours,  and  finally 
a  return  to  normoxia.  The  model  predicts  a  decrease  in  brain  oxygen  levels,  a  corresponding 
increase  in  HIF,  and  a  resulting  increase  in  transcription  factor. 

Simulations  such  as  these  can  be  used  to  describe  and  interpret  experimental  data  in  terms  of 
underlying  mechanisms.  Such  experiments,  if  quantitative,  can  also  be  used  to  calibrate  and 
further  parameterize  the  model. 
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Figure  8:  Simulated  Kinetics  of  O2  in  Arterial  Blood,  Brain  Blood,  and  Brain  Tissue. 

Normoxia  (21  percent),  hypoxia  (7  percent),  and  return  to  normoxia.  Oxygen  concentrations  are 
normalized  to  levels  under  normoxia  (1.0  for  21  percent  O2  in  inhaled  air),  assuming  a  decrease 
in  hemoglobin  binding  (and  arterial  concentration)  at  7  percent  O2  by  about  37  percent  (see 
Figure  4A). 
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Figure  9.  Simulated  Kinetics  of  HIF-la  in  Brain  Cytosol  and  Nucleus,  Together  with 
Transcription  Signal.  Normoxia  (21  percent),  hypoxia  (7  percent),  and  return  to  normoxia. 


4.2  Dose- Response  Simulations 

Although  the  model  can  be  extended  to  include  gene  transcription  (see  Section  5.1  below),  we 
can  analyze  steady-state  (SS)  gene  transcription  products  without  a  detailed  time-course  model. 
For  example,  Schoch  et  al.  (2002)  measured  the  amount  of  VEGF  protein  and  mRNA  in  the 
brain  of  mice  following  normobaric  hypoxia  at  12  percent  -  6  percent  oxygen  for  24  hours  in 
vivo  using  an  ELISA  (enzyme-linked  immunosorbent  assay)  specific  for  murine  VEGF  and 
northern  blots,  respectively.  These  data  are  shown  as  the  red  and  blue  plots  in  Figure  10. 

The  mRNA  and  protein  levels  in  the  SS  will  be  given  by  the  ratio  of  their  production  and 
degradation  rates.  It  is  reasonable  to  assume  that  degradation  rates  are  relatively  constant,  and  in 
particular  are  not  dependent  on  the  degree  of  hypoxia  or  hyperoxia.  Production  rates,  however, 
may  reasonably  be  assumed  to  be  proportional  to  the  degree  of  hypoxia  (linear  dose/response),  at 
least  within  reasonable  limits.  The  SS  levels  of  mRNA  and  protein  will  thus  also  be  proportional 
to  the  (hypoxic)  dose.  If  we  further  assume  that  in  the  SS,  levels  of  protein  and  mRNA  are 
proportional  to  the  transcription  signal  given  from  our  model  TC  (double  the  transcription  signal 
leads  to  doubling  of  the  SS  levels  of  protein  and  message),  we  may  fit  the  Schoch  et  al.  (2002) 
dose/response  data  as  shown  in  Figure  8,  where  we  have  plotted  the  ratio  of  protein  and  message 
at  12  percent,  10  percent  8  percent  and  6  percent  O2,  relative  to  normoxia  (20  percent  O2), 
against  the  respective  O2  concentrations.  Our  results  indicate  that  our  model  is  able  to  reasonably 
fit  the  changing  ratios  of  both  VEGF  protein  and  mRNA  as  a  function  of  hypoxia,  although  the 
fitted  value  of  prolyl  hydroxylase  enzyme  activity  (PHD)  was  lower  by  a  factor  of  about  1 .7/4  ~ 
0.4  for  the  mRNA  fit  (Figure  10). 
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Figure  10.  Transcription  Signal  (TC)  Fits  to  mRNA  and  VEGF  Data.  The  Y  axis  shows  the 
ratio  of  protein  and  message  at  12  percent,  10  percent  8  percent  and  6  percent  O2,  relative  to 
normoxia  (20  percent  O2),  against  the  respective  O2  concentrations.  Red  curve  represents 
measured  VEGF  protein  levels,  blue  curve  mRNA  levels  for  VEGF.  Gray  lines  represent  model 
predictions,  for  different  values  of  prolyl  hydroxylase  activity  (PHD)  optimized  to  fit  each  data 
set.  Error  bars  represent  standard  deviations  for  the  ratios.  Data  from  Schoch  et  al.  (2002). 


5.0  CONCLUSIONS  AND  FUTURE  WORK 


5.1  Gene  Transcription  and  Longer  Term  Effects 

The  model  can  be  extended  to  include  a  gene  transcription  signal  and  gene  transcription  itself  as 
output  from  the  HIF- la  pathway.  In  such  an  extended  model,  the  output  from  the  current  model 
(HIF-la  binding  to  the  HRE  and  the  resulting  transcription  signal)  would  provide  an  input  to  a 
model  of  gene  expression,  such  as  one  provided  by  Chen  et  al.,  1999  or  Ay  and  Amosti,  2011 
(see  Figure  11),  thereby  linking  to  angiogenesis  and  finally  to  alterations  in  capillary  density 
(LaManna  et  al.,  2004). 
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Feedback  Control 


n  The  number  of  genes  in  the  genome; 

r  mRN  A  concentrations,  ^-dimensional  vector- valued  functions  of  f; 

p  Protein  concentrations,  n-dimensional  vector-valued  functions  of  t; 

/( p)  Transcription  functions,  n-dimensional  vector  polynomials  on  p; 

L  Translational  constants,  n  x  n  non-degenerate  diagonal  matrix; 

V  Degradation  rates  of  mRNAs;  n  X  n  non-degenerate  diagonal  matrix; 

U  Degradation  rates  of  Proteins,  n  x  n  non-degenerate  diagonal  matrix; 

HRE  Hypoxia  Responsive  Element 

Figure  11.  Schematic  of  a  Simple  Differential  Equation  Based  Model  of  Gene  Expression. 

Based  on  Chen  et  al.  (1999). 


5.2  Capillary  Densities  and  Angiogenesis 

VEGF  and  other  proteins,  released  from  cells  as  a  result  of  local  hypoxic  conditions,  stimulate 
chemotaxis  and  proliferation  in  endothelial  cells  during  capillary  sprouting.  Experimental  data 
are  available  that  explore  the  time-course  of  angiogenesis  and  changes  in  capillary  density 
following  extended  periods  of  hypoxia,  and  return  to  normoxia.  Figure  12  (A)  shows  HIF-la 
measured  in  the  cortex  of  rats  exposed  to  hypoxic  conditions  for  6  hours  to  21  days,  relative  to 
normoxic  controls  (0.5  atmospheres,  equivalent  to  10  percent  normobaric  oxygen)  (Chavez  et  al., 
2000).  Figure  12  (B)  shows  capillary  density  in  the  rat  cerebral  cortex  increased  by  60  percent 
after  3  weeks  of  hypoxia  (0.5  atmospheres,  equivalent  to  10  percent  normobaric  oxygen)  and  that 
it  progressively  decreased  to  prehypoxic  values  after  3  weeks  of  normoxic  recovery  (Pichiule 
and  Fa  Manna  2002).  The  latter  is  likely  due  to  apoptosis  of  cells  associated  with  excess 
capillary  density  that  is  no  longer  required  under  normoxic  conditions.  Both  capillary  density 
increases  (due  to  angiogenesis)  and  decreases  (due  to  apoptosis)  have  similar  time  scales,  taking 
days  or  weeks  to  adjust,  making  them  unlikely  candidates  to  mediate  more  rapid  changes  in 
hypoxia  or  periods  of  rapid  intermittent  hypoxia. 
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Figure  12.  HIF-1  a  and  Capillary  Density  in  Hypoxic  Rat  Cerebral  Cortex.  (A)  HIF-la 
Measured  in  the  Cortex  of  Rats  Exposed  to  Hypoxic  Conditions:  Exposure  is  for  6  hours  to  21 
days,  relative  to  normoxic  controls  (0.5  atmospheres,  equivalent  to  10  percent  normobaric 
oxygen)  (Chavez  et  al.,  2000).  (B)  Capillary  Density  in  the  Rat  Cerebral  Cortex:  This  increased 
by  60  percent  after  3  weeks  of  hypoxia  (0.5  atmospheres,  equivalent  to  10  percent  normobaric 
oxygen)  and  that  it  progressively  decreased  to  prehypoxic  values  after  3  weeks  of  normoxic 
recovery  (adapted  from  Pichiule  and  La  Manna  2002,  based  on  photomicrographs  of  GLUT-1- 
stained  sections). 


Models  of  angiogenesis  have  been  developed  (see  Qutub  et  al.,  2009  for  an  overview).  Basically, 
the  process  is  as  follows  (Qutub  et  al.,  2009): 

1)  HIF-la  is  upregulated  in  a  hypoxic  cell.  HIF-la  activates  the  transcription  of  VEGF, 
which  is  then  secreted  by  the  cell. 

2)  VEGF  -  VEGF-receptor  binding  on  the  capillary  surface. 

3)  Vessel  permeability  changes. 

4)  An  activated  endothelial  cell  (the  tip  cell)  starts  to  break  down  the  basement  membrane. 

5)  Stalk  cells  proliferate  behind  the  tip  cell. 

6)  The  leading  edge  of  the  moving  sprout  releases  matrix  metalloproteinases,  which 
proteolyze  the  surrounding  extracellular  matrix,  allowing  the  cell  to  migrate. 

Such  a  processes  is  necessarily  multi- scale  in  nature  -  and  multiscale  modeling  is  an  essential 
tool  to  truly  understand  a  process  of  such  biological  complexity. 


5.3.  Inflammation 

HIF  can  have  both  pro-  and  anti-inflammatory  effects,  and  cross-talk  between  HIF  and 
inflammatory  pathways  may  be  critical  for  understanding  white  matter  hyperintensities  (WMH) 
in  U-2  pilots,  since  the  latter  are  likely  associated  with  neuroinflammation  (Cavadas  et  al.,  2013; 
McGuire  et  al.,  2013).  VEGF  induces  angiogenesis,  increased  vascular  permeability,  and 
inflammation.  VEGF  has  been  shown  to  induce  angiogenesis  when  infused  continuously  into 
adult  rat  brain  tissue;  in  addition,  VEGF  has  been  shown  to  enhance  permeability  in  brain 
vasculature  (Vezzani,  2005). 
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Pilots  of  the  high-altitude  U-2  reconnaissance  aircraft  undergo  repeated  exposure  to  hypobaric 
conditions,  and  as  a  result  may  suffer  brain  lesions  (white  matter  hyperintensities)  and  cognitive 
deficits  (McGuire  et  al.,  2014).  All  of  the  WMH  observed  in  pilots  were  regionally  uniformly 
distributed  and  located  in  deep  white  matter  rather  than  in  the  cerebral  cortex  (McGuire  et  al., 

2013) ,  while  there  is  a  significant  regional  heterogeneity  in  the  distribution  of  subcortical  WMH 
in  normal  aging.  The  majority  (60  percent-80  percent)  of  WMH  are  found  in  the  frontal  area, 
presumably  because  its  high  metabolic  demand  makes  it  more  vulnerable  to  age-related 
cerebrovascular  disorders.  A  more  uniform  distribution  of  WMH  is  a  hallmark  finding  in  many 
neuroinflammatory  disorders  and  traumatic  brain  injury  and  may  be  used  to  gauge  disease 
severity  and  progression,  suggesting  that  hypobaric  exposure  produces  white  matter  damage 
different  from  that  occurring  in  normal  aging.  Other  groups  with  increased  WMH  burden  include 
high-altitude  mountain  climbers,  in  which  WMH’s  are  attributed  to  a  combination  of  hypoxia 
and  hypobaria  (McGuire  et  al.,  2013),  and  high  altitude  chamber  operators  (McGuire  et  al., 
2014b).  The  latter  study  provides  strong  evidence  that  non-hypoxic  hypobaric  exposure  may 
induce  subcortical  WMHs  in  a  young,  healthy  population  lacking  other  risk  factors  for  WMHs 
(McGuire  et  al.,  2014b). 

The  circulation  of  U-2  pilots  (and  high  altitude  chamber  operators)  is  purged  of  nitrogen 
(degassed)  prior  to  each  mission  with  100  percent  oxygen  for  1  hour  (McGuire  et  al.,  2014).  This 
may  not,  however  be  sufficient,  since  gases  readily  diffuse  across  the  mucosa.  The  direction  of 
diffusion  is  dictated  by  the  partial  pressure  of  that  gas  in  blood  versus  luminal  contents.  For 
methane  and  hydrogen,  diffusion  is  always  out  of  the  lumen  into  blood.  Nitrogen  and  CO2 
diffuse  in  either  direction,  depending  on  specific  conditions  within  the  individual.  We  thus 
postulate  that  microinfarcts  can  be  produced  by  nitrogen  bubbles,  inducing  microvascular 
occlusions  (McGuire  et  al.,  2013).  Damage  produced  by  interaction  between  microemboli  and 
cerebral  tissue,  may  lead  to  thrombosis,  coagulation,  inflammation,  and/or  activation  of  innate 
immune  response.  Injury  may  thus  be  produced  by  microemboli  entering  cerebral  circulation  at 
random. 

According  to  McGuire  et  al.,  (2013),  three  potential  sources  of  microemboli  should  be 
considered:  microbubbles  of  gas,  presumably  nitrogen;  platelet-based  thrombi;  and 
microparticles. 

•  Microbubbles:  symptoms  of  decompression  sickness  classically  are  believed  secondary  to 
nitrogen  gas  bubbles  exerting  direct  pressure  on  tissues,  blocking  small  arteriolar  vessels, 
and  interacting  with  blood  proteins 

•  Microthrombi:  occlusion  of  small  cerebral  vessels  by  platelet  thrombi  produced  by 
accelerated  coagulation  of  blood  in  the  presence  of  venous  nitrogen  gas  bubbles 

•  Platelet  aggregation. 

Thus,  vascular  damage  or  tissue  damage  caused  by  these  processes  induces  (and  is  exacerbated 
by)  an  inflammatory  response.  In  addition,  IL-ip  (and  potentially  other  cytokines)  released  by 
microglia  and  macrophages  impacts  long  term  potentiation  (LTP)  and  memory  formation 
(Lynch,  2010),  potentially  leading  to  the  observed  cognitive  issues  in  U-2  pilots  (McGuire  et  al., 

2014)  (see  Figure  13). 
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Furthermore,  and  more  generally,  periods  of  hypoxia,  followed  by  normoxia  (or  hyperoxia) 
resulting  from  various  environmental  factors  that  may  be  encountered  by  pilots  (such  as  rapid 
acceleration  or  oxygen  supply  issues),  may  lead  to  brain  tissue  damage  due  to  oxidative  stress 
(and  inflammatory)  mechanisms  associated  with  reperfusion  injury.  In  such  cases,  when  the 
tissue  blood  supply  returns  to  normal  after  a  period  of  ischemia  or  lack  of  oxygen,  the  prior 
hypoxic  conditions  creates  a  condition  in  which  the  restoration  of  circulation  (reperfusion) 
results  in  inflammation  and  oxidative  damage  through  the  induction  of  oxidative  stress, 
potentially  leading  to  cognitive  issues  (Carden  &  Granger,  2000). 


Figure  13.  Relation  between  Brain  Tissue  Hypoxia,  Neuroinflammation  and  Performance 
Degradation 


Cytokine  release  following  hypoxia  has  been  observed  experimentally.  For  example,  Kozak  et 
al.  (2006)  demonstrated  that  mice  exposed  for  7  days  to  hypoxia  (11  percent  O2)  resulted  in 
significant  elevation  of  plasma  IL-6  (see  Figure  14),  inducing  sickness  behavior,  a  set  of  adaptive 
responses  of  the  host  to  severe  infections  and  inflammation.  It  includes,  among  others,  the 
thermoregulatory  responses  such  as  regulated  increase  (fever)  and/or  decrease  of  body 
temperature  (Tb),  decrease  of  motor  activity  (lethargy),  and  loss  of  appetite  (Kozak  et  al.,  2006). 
Such  time-course  data  may  ultimately  be  described  by  a  model  that  includes  HIF-la  -  induced 
gene  transcription  and  further  pro-inflammatory  processes. 
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□  Normoxia  (n=6/time  point) 


Time  (h) 


Figure  14.  Changes  of  Plasma  IL-6  Levels  over  Time.  Swiss  Webster  mice  were  exposed  to 
hypoxia  (11  percent  O2;  closed  bars)  and  normal  air  (open  bars).  Plasma  IL-6  levels  are 
expressed  in  pg/mL  and  time  in  hours  (h).  From  Kozak  et  al.  (2006). 


5.4  Carotid  Body  and  Other  Issues 

Application  of  the  HIF-la  model  to  the  carotid  body  is  more  speculative  (see  Introduction). 
However,  “hydroxylation  of  HIF  in  the  presence  of  O2  occurs  in  a  few  minutes,  hence  it  is 
conceivable  that  02-dependent  hydroxylases  could  also  modulate  ion  channels  and  thus 
participate  in  the  acute  responses  to  hypoxia”  (Lopez-Barneo  el  al.,  2008).  The  carotid  body 
response  on  heart  rate,  breathing  rate,  clearly  needs  to  be  rapid.  Rapid  response  pathways 
potentially  mediated,  at  least  in  part,  by  the  HIF-la  pathway,  include  nitric  oxide  synthase 
(NOS)  activation  (Jung  et  al.,  2000),  vasodilation,  glucose  transport,  and  metabolism.  If  we 
assume  that  oxygen  sensing  is  again  mediated  by  the  intracellular  levels  of  HIF-la,  then  we  have 
the  scheme  outlined  in  Figure  15. 

Transport  of  HIF-la  to  nucleus  takes  place  via  microtubules  -  filamentous  intracellular 
structures  that  are  responsible  for  various  kinds  of  movements  in  all  eukaryotic  cells. 
Microtubules  are  involved  in  nucleic  and  cell  division  (chromosome  separation  in  mitosis), 
organization  of  intracellular  structure,  and  intracellular  transport.  Carbonaro  et  al.  (2012)  showed 
that  HIF-la  protein  associates  with  polymerized  microtubules  and  is  transported  to  the  nucleus 
with  the  aid  of  the  dynein  motor  protein.  Microtubule-targeting  drug  (MTD)  treatment  impaired 
HIF-la  protein  nuclear  translocation,  which  significantly  down-regulated  HIF  transcriptional 
activity.  Understanding  of  the  microtubule-dependent  HIF-la  regulation  is  essential  given  the 
importance  of  HIF-la  in  tumor  biology,  particularly  tumor  angiogenesis,  and  the  widespread  use 
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of  MTDs  in  clinical  oncology  (Carbonaro  et  al.,  2012).  Because  of  the  importance  of 
microtubules  in  a  myriad  of  other  transport  processes,  particularly  regarding  transport  of 
receptors  and  ion  channels  to  membranes  in  cognition  and  memory  formation,  disruption 
occurring  with  MTDs  (and  other  stressors)  clearly  has  far-reaching  implications.  The 
involvement  of  microtubules  in  both  hypoxia  signaling  and  cognition  also  suggests  the  potential 
for  significant  cross-talk  between  these  processes. 


Figure  15.  Schematic  of  a  Potential  Model  for  Oxygen  Sensing  in  the  Carotid  Body 


In  many  cases  oxygen  delivery  issues  are  concurrent  with  ambient  diverse  chemical  exposure 
issues.  For  example,  toluene  is  a  known  contamination  component  that  comes  through  the 
aircraft’s  on-board  oxygen  generation  system,  and  may  potentially  impact  oxygen  delivery 
(Martin  et  al.,  2012).  Previous  in-house  experiments  have  investigated  the  effect  of  altitude  on 
the  distribution  of  toluene  in  rats  by  providing  a  better  understanding  of  how  these  stressors 
affect  the  brain  at  the  cellular  and  molecular  level.  Ongoing  studies  with  toluene  combined  with 
hypoxia  are  investigating  the  effects  of  altitude,  chemical  exposure  and  altered  percentage  of 
oxygen  on  the  function  of  neural  cells  in  the  brain  by  measuring  changes  in  the  electrical  patterns 
of  brain  cells.  Toluene  also  may  potentially  interact  with  hemoglobin,  altering  its  binding 
characteristics  and  thereby  the  oxygen  carrying  capacity  of  blood  (Chushak  et  al.,  2015). 

Oxygen  delivery  to  tissues  has  been  significantly  studied;  the  current  oxygen  transport  model  is 
necessarily  a  simplification  of  a  complex  process.  Improvements  to  the  oxygen  transport  model, 

24 


Distribution  A:  Approved  for  public  release;  distribution  unlimited  (PA  Case  No  88ABW-2017-6197,  11  Dec  2017) 


based  on  this  extensive  literature,  could  include  more  dynamic  detail,  such  as  a  4- stage  binding 
approach  with  Kon  and  K0ff  rates.  Delivery  to  brain  also  depends  critically  on  local  cerebral  blood 
flows,  which  respond  to  metabolic  demands  based  on  brain  functional  activity.  Control  of 
cerebral  blood  flow  is  complex  and  only  partially  understood.  Relevant  processes  include 
cerebral  pressure  autoregulation,  which  maintains  a  constant  flow  despite  changing  cerebral 
perfusion  pressure.  The  brain  is  also  able  to  vary  blood  flow  to  match  metabolic  activity  (flow- 
metabolism  coupling).  Finally,  a  network  of  perivascular  nerves  also  serves  to  regulate  cerebral 
blood  flow  (neurogenic  regulation)  (Peterson  et  al.,  2011). 

We  have  treated  PHD  in  the  model  as  a  single  enzyme  -  in  reality,  3  PHD  isoforms  have  been 
identified,  that  contribute  to  the  regulation  of  HIF-la  (Appelhoff  et  al.,  2004).  It  is  hypothesized 
that  their  differing  cell  specific  behaviors  may  lead  to  flexibility  in  the  HIF-la  response.  Once 
these  behaviors  are  more  fully  characterized,  a  more  subtle  response  to  hypoxia  could  be 
incorporated  into  the  model. 
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APPENDIX  A.  MODEL  CODE:  COMBINED  OXYGEN  TRANSPORT  AND  HIF-1A 

MODEL 


Note  that  the  HIF-la  model  was  originally  encoded  in  Berkeley  Madonna  software  and 
presented  in  a  previous  Technical  Report  (Robinson  et  al.,  2016) 

#  02  transport  and  HIF  response  model 

#  Version  7  .  1 

#  8  September  2017 

#  C.  Eric  Hack  and  Peter  J.  Robinson 

#  Kinetic  model  of  oxygen  delivery  to  the  brain  and  HIFla  signaling. 

#  lines  beginning  with  '#'  are  comments,  and  not  executed 

#  clean  the  workspace 
rm ( list=ls ( ) ) 

#  load  the  solver 
library (deSolve) 


###################################################################### 

#  Constants 

#  - Initial  conditions  for  the  amounts 

ATBO  <-  4 

ATTO  <-  43 
ABrBO  <-  0.01 
ABrTO  <-  1 
HIF0  <-  le-9 
HIFN0  <-  le-14 
TC0  <-  5 . 3e- 1 1 


# - Chemical  parameters 

MW  <-32.0  #  Molecular  weight  (g/mole)  of  02 

MW2  <-44.0  #  Molecular  weight  (g/mole)  of  metabolite  (C02) 


# - Physiological  parameters 


QPC  <-  14. 
QCC  <-  14. 
QBrC  <-  0.10 
1983) 

QTC  <-  0.9 


#  Alveolar  ventilation  rate  (L/hr) 

#  Cardiac  output  (L/hr) 

#  Fractional  blood  flow  to  brain  (L/hr;  Gerlowski  and  Jain, 

#  Fractional  blood  flow  to  rest  of  body  (L/hr) 


BW  <- 
VBrC  <- 
VBrBC  <- 
VTC  <- 
perfused) 
VBld  <- 
VTBC  <- 


70 

0 . 02 
0 . 05 
0 .89 

0 . 075 
VBld  - 


#  Body  weight  (kg) 

#  Fraction  brain  tissue  (kg  or  L) 

#  Fraction  of  brain  that  is  blood 

#  Fraction  rest  of  body  tissue  (kg  or  L;  total  ~91% 
#  Fraction  blood  in  while  body 

VBrC*VBrBC  #  Fraction  of  rest  of  body  that  is  blood 


#  Partition  coefficients 
PB  <-  0.03  #  Blood:air 

PBr  <-  1  #  Brain:blood 

PT  <-  1  #  Rest  of  body  tissues : blood 

PABr  <-0.01  #  Permeability  X  Area  cross 

blood  and  brain  tissue 

PAT  <-1.01  #  Permeability  X  Area  cross 

blood  and  tissue  cells 


product  between  brain 
product  between  tissue 
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#  Metabolism  parameters 

VMAXBrC  <-  1  #  Brain  Maximum  rate  of  metabolism  (mg/hr) 

VMAXTC  <-  1  #  Tissue  Maximum  rate  of  metabolism  (mg/hr) 

KM  <-  1  #  Half-max  concentration  (mg/L) ;  not  scaled  by  BW 

#  - Inhalation  exposure  parameters 

P  <-  1  #  Pressure  (atm) 

PmmHg  <-  P*760  #  Pressure  (mmHg) 

F02  <-  0.21  #  Fraction  02  in  outside  air 

PH20  <-  47  #  Water  vapor  pressure  in  alveolus  (mmHg) 

humid  <-  (PmmHg  -  PH2O)/760  #  Note  that  this  is  about  6%  reduction  from  760 
mmHg,  and  agrees  with  Karius . 

#  (Karius  DR,  "Respiratory  Adaptations  in  Health  and 
Disease:  Calculating  the  Alveolar 

#  Oxygen  &  the  A-a  02  gradient",  a  Kansas  City  University 
of  Medicine  and  Biosciences  online  lecture) . 

CC  <-  FALSE  #  Default  to  open  chamber 
NRATS  <-  3 .  #  Number  of  rats  (for  closed  chamber) 

KLC  <-  0 .  #  First  order  loss  from  closed  chamber  (/hr) 

VCHC  <-9.1  #  Volume  of  closed  chamber  (L) 

TCHNG  <-12.0  #  Length  of  inhalation  exposure  (hrs) 

#  - Timing  commands 

#?  How  to  incorporate  multiple  changes  in  exposure  in  R? 

#?  Are  the  good,  built-in  forcing  functions,  like  step  functions? 

TSTOP  <-24.0  #  Length  of  experiment  (hrs) 

POINTS  <-  100  #  Number  of  points  in  plot 

CINT  <-  TSTOP/POINTS  #  Communication  interval  (for  output,  hrs) 


#  - Uptake  of  inhaled  02  based  on  Hgb  saturation  curve 

#  The  hbg  binding  parameters  are  set  in  the  Fit.CB.R  script. 

#  Retained  here  for  documentation  purposes. 


Bmax  <-  1  #  fraction  Hgb  occupied  by  02 

KD  <-  42  #  mmHg 

n  <-  2.4  #  Hill  coefficient  for  02  binding  to  Hgb 

hgb  <-  150  #  g  Hgb  /  L  blood 

#  - Ideal  gas  law  parameters  for  conversion  between  mL  02  and  mg  02  (CA) 

#  The  conditions  are  for  STPD,  since  the  mL  represents  what  would  be  the 
volume  at  STPD. 

R  <-  0.082056  #  Universal  gas  constant  (L*atm/K/mol) 

Temp  <-  273.15  #  Standard  temperature  (K) 

BTemp  <-  Temp  +37  #  Body  temperature  (K) 


#  - Conversion  factor  for  STPD  to  BTPS 

StoB  <-  (760)/ (PmmHg  -  PH20)  *  (BTemp)/ (Temp)  #  =  863/ (PmmHg  -  47) 

#  - Dummy  parameter  for  CA  and  CV  extraction 


kbig  <-  1 . 0e9 

#  Get  relationship  between  CF  (free)  and  CB  (bound) . 
source ("Fit.CB.R") 

#  The  parameters  are  assigned  by  the  script  file. 


#  Bmax  <-  coef(fit)  [1] 

#  KD  <-  coef(fit)  [2] 

#  n  <-  coef  (fit)  [3] 

#  - HIF  submodel-specific  parameters 

#  most  parameters  are  from  Nguyen  et  al .  (2013) 

Kphd  <-  162  #  /hr 

Tmax  <-  le-3/le3*3600  #  uM/hr 

KT  <-  1  #  uM 

Khif  <-  1.8e-2  #  uM/hr 

Kbm  <-32.4  #  /hr 
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#  Kbm*HRE 


32 . 4/hr 


HRE  <-l 

KdegN  <-  1 

# 

/hr 

Ktc  <-  0.36 

# 

/hr 

PHD  <-  5000 

# 

uM 

KmHIF  <-  0.1 

Km02  <-  250 

# 

uM 

#  - Modify  the  inhaled  02  fractions  and  duration  (change  times  are 

calculated) 

#  02  fractions 

Oxy  <-  c  (0 . 21,  0.105,  0.21,  0.21) 

#  Change  times  and  durations 

Dur  <-  c(l,  6,  1,  Inf) 

tau  <-  c(0,  Dur[l]  +  Dur[2],  sum ( Dur [ 1 : 3 ] ) ) 

###################################################################### 

#  Modify  constants 

#  Source  a  file  with  modifications  to  the  default  constant  parameter  values. 

#  This  needs  to  be  done  before  calculations  based  on  the  constants. 

#  For  example,  load  a  rat  or  human  physiological  parameter  file, 

#  or  change  the  dosing  parameters,  or  change  the  chemical. 

#  Can  have  multiple  R  scripts  called  here. 

#BUT  BE  CAREFUL  OF  THE  ORDER! 

#(that  is,  when  changing  the  same  parameter  in  multiple  files) 


#  Set  species-specific  parameters 

source ("rat.R")  #  sets  BW  to  0.3,  VFC  to  0.09 

#source ( "human . R" ) 


#  Set  chemical-specific  parameters 

#source ( "Styrene . R" )  #  placeholder  for  future  addition  of  other  chemicals 


#  Set  exposure  scenario 
source ( "hypoxia2 . R" )  # 
#source ( "modify . R" )  # 
#source ("johnson.R")  # 


load  a  particular  hypoxic  event 

placeholder  for  other  parameter  modifications 

simulate  the  Johnson  study 


##################################################################### 

#  Scaled  or  Adjusted  or  Calculated  parameters 


#  convert  from  fraction  02  in  atmosphere,  and  decrease  for  humidification. 
CONC  <-  P*F02*1 . 0e6*humid  #  Inhaled  concentration  (ppm) 


if  (CC)  {  #  Closed  chamber  simulation 
RATS  <-  NRATS 
KL  <-  KLC 
}  else  { 

RATS  <-  0. 

KL  <-  0 .  #  (Turn  off  chamber  losses  so  concentration  remains  constant) 


VCH  <-  VCHC-RATS*BW  #  Net  chamber  volume  (L) 

AI0  <-  CONC*VCH*MW/24450 .  #  Initial  amount  in  chamber  (mg) 


#  - Allometric  scaling 

#  Flows  (blood  and  respiration) 

QC  <-  QCC*BW* *0.74#  Cardiac  output  (L/hr) 

QP  <-  QPC*BW**0.74  #  Alveolar  ventilation  (L/hr) 
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QBr  <-  QBrC*QC  #  Brain  blood  flow  (L/hr) 

QT  <-  QTC*QC  #  Rest  of  tissues  blood  flow  (L/hr) 

#  Tissue  Volumes 
VBr  <-  VBrC*BW 
VBrB  <-  VBrBC*VBr  # 

VBrT  <-  VBr  -  VBrB 
VT  <-  VTC*BW 
VTB  <-  VTBC*VT 
VTT  <-  VT  -  VTB 

#  Metabolism 

VMAXBr  <-  VMAXBrC*BW* *  0 . 7  #  Maximal  rate  (mg/hr) 

VMAXT  <-  VMAXTC*BW**0 . 7  #  Maximal  rate  (mg/hr) 

#  Conversion  factors 

nMps2mgph  <-  VBrT*3600*MW/le6  #  convert  nM/s  to  mg/hr 

#  nMps2mgph  needed  to  use  the  HIF-PHD-02  equation  in  the  brain  tissue  rate  eqn. 


###################################################################### 

#  Define  the  model/system  of  ODEs 

dynamic  <-  function  (times,  state,  parms)  { 
with (as . list (c (state,  parms)),  { 

#  - DYNAMIC:  Things  that  change  with  time  go  in  this  function  definition 

#  - Compartments 

#  Brain  tissue 

#  Brain  blood 

#  Alveolar  gas  exchange  (Cl,  CA,  CX) 

#  Arterial  blood  Hgb  (CAF,  CAB) 

#  Rest  of  body  tissue  (cells) 

#  Rest  of  body  blood 

#  - Inhalation 

#  idx  is  exposure  index  variable 

idx  <-  0 

len  <-  length (tau) 

for  (i  in  l:len)  {  idx  <-  idx  +  as . numeric ( times  >=  tau[i])  } 
frax  <-  Oxy [ idx] /Oxy [ 1 ] 

Cl  ZONE  <-  as . numeric (CC  |  (times  <  TCHNG) )  #  0  if  CC=False  and  t>=TCHNG,  1 
otherwise 

Cl  <-  AI/VCH*CIZONE*frax  #  Concentration  in  alveolar  air 

(mg/L) 

CP  <-  CI*24450./MW  #  Alveolar  Air  Concentration  (ppm) 

#  - Body,  or  Tissues  (T) 

CTT  <-  ATT/VTT  #  tissue/cells 
CVT  <-  ATB/VTB  #  tissue  blood 

#  - Brain 

CBrT  <-  ABrT/VBrT 
CVBr  <-  ABrB/VBrB 

#  - CV  =  Mixed  venous  blood  concentration  (mg/L) 

CV  <-  (QT*CVT  +  QBr*CVBr ) / QC 

#  - PA02  =  alveolar  P02  (mmHg) 

PA02  <-  P* 7 60*Oxy [ idx] *humid  #  alveolar  02  (mmHg) 

#  - CA  =  Concentration  in  arterial  blood  (mg/L) 
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#  Brain  total  (L) 

Brain  blood  (L) 

#  Brain  tissue  (L) 

#  Rest  of  tissues  (L) 

#  Tissue  blood  (L) 

#  Tissue  tissue/cells  (L) 


#  Implicit  solve  for  CAF .  Find  the  zero  of  the  function  using  uniroot () . 

f  <-  function (x)  {QP*(CI  -  x/PB)  +  QC* (CV  -  Bmax*x**n/ (KD**n  +  x**n)  -  x) } 

soln  <-  uniroot (f,  c(0,300)) 

CAF  <-  soln$root 

#  Now  compute  bound  02 

CAB  <-  Bmax*CAF**n/ (KD**n  +  CAF**n) 

#  Total  arterial  02 
CA  <-  CAF  +  CAB 

#  Exhaled  02 

CX  <-  CAF/PB  #  Concentration  in  exhaled  air  (mg/L) 

CXPPM  <-  ( 0 . 7 *CX+0 . 3*CI ) *2 4 4 50 . /MW  #  Concentration  (ppm),  30%  dead  space 

assumed 

#  - A  =  dummy  for  CA 

dA  <-  kbig* (CA  -  A) 

#  - V  =  dummy  for  CV 

dV  <-  kbig* (CV  -  V) 

#  - AI  =  Amount  in  inhaled  air  (mg/L) 

dAI  <-  RATS*QP* (CAF/PB-CI)  -  (KL*AI)  #  Rate  of  change  in  air  amount 
(mg/hr) 

dAInh  <-  QP*CI  #  Rate  inhaled  (mg/hr) 

#  - AX  =  Amount  exhaled  (mg) 

dAX  <-  QP*CX  #  Rate  exhaled  (mg/hr) 

#  - ABr  =  Amount  in  brain  blood  (mg) 

brup  <-  PABr* (CVBr  -  CBrT/PBr)  #  oxygen  diffusion  from  brain  blood  to 
tissue  (mg/hr) 

dABrB  <-  QBr* (CA  -  CVBr)  -  brup 

#  - AT  =  Amount  tissues  in  rest  of  body  (mg) 

dATB  <-  QT* (CA  -  CVT )  -  PAT* (CVT  -  CTT/PT)  #  tissue  blood 

dATT  <-  PAT* (CVT  -  CTT/PT)  -  VMAXT*CTT/ (KM  +  CTT ) 

#  - AM  =  Amount  metabolized  in  brain  and  rest  of  body  (mg) 

brmet  <-  VMAXBr*CBrT/ (KM  +  CBrT ) 
tmet  <-  VMAXT*CTT/ (KM  +  CTT) 

hifmet  <-  Kphd*PHD*HIF*CBrT/ ( (HIF+KmHIF) * (CBrT+Km02) ) 
dAM  <-  brmet  +  tmet  +  hifmet  #  total  metabolism 


# - HIFla  submodel 

#  Oxygen  diffuses  in  from  blood  and  is  removed  from  cytosol  via 
respiration  and  hydroxylation  of  HIF 

dABrT  <-  brup  -  brmet  -  hifmet*nMps2mgph  #  mg/hr 

#  HIF  is  synthesized,  and  removed  from  cytosol  by  hydroxylation  and 
transport  to  nucleus 

hiftrans  <-  Tmax*HIF/ (HIF+KT)  #  uM/hr 

dHIF  <-  Khif  -  hifmet  -  hiftrans  #  uM/hr 

#  HIF  is  removed  from  nucleus  by  degradation  and  binding  to  hypoxia 
responsive  element  (HRE) 

dHIFN  <-  hiftrans  -  KdegN*HIFN  -  Kbm*HIFN*HRE  #  uM/hr 

#  Transcription  signal  (TC)  creation  by  binding  of  HIF  to  HRE,  and 
degradation 

dTC  <-  Kbm*  H I FN  *  HRE  -  Ktc*TC  #  uM/hr 
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# - model  function  must  return  a  list  for  the  ode  solver 

return (list (c (dAI,  dAInh,  dAX,  dABrB ,  dATB ,  dATT,  dAM,  dABrT ,  dHIF,  dHIFN, 
dTC,  dA,  dV) ) ) 

}) 

} 


##################################################################### 

#  Integration 

#  Parameters  for  integrator 

#  Using  global  parameters  for  now, 

#  so  parms  is  a  dummy  to  satisfy  the  input  requirements  of  our  derivative 
function . 

parms  <-  1 

#  States 

state  <-  c (AI  =  AIO, 

Alnh  =  0, 

AX  =  0, 

ABrB  =  ABrBO, 

ATB  =  ATBO , 

ATT  =  ATTO, 

AM  =  0, 

ABrT  =  ABrTO, 

HIF  =  HIF0 , 

HIFN  =  HIFNO , 

TC  =  TCO , 

A  =  200, 

V  =  200) 

#  Integration  times 

times  <-  seq ( 0 ,  TSTOP,  by  =  CINT) 

#  Call  the  integrator 

#out  <-  ode(y  =  state,  times  =  times,  func  =  dynamic,  parms  =  parms, 
atol=l . Oe-12 ,  rtol=1.0e-9) 

out  <-  ode(y  =  state,  times  =  times,  func  =  dynamic,  parms  =  parms) 
y  <-  as . data . frame (out) 


##################################################################### 

#  - Mass  balance  equations 

#  - TMASS  =  total  mass  stored  (mg) 

ABr  <-  out [, "ABrT"]  +  out [, "ABrB"] 

AT  <-  out [, "ATT"]  +  out [, "ATB"] 

TMASS  <-  ABr  +  AT 

T MAS SO  <-  ABrBO  +  ABrTO  +  ATBO  +  ATTO 

#  - OUT  =  amount  removed  by  metabolism  and  exhalation  (mg) 

OUT  <-  out [ , "AM" ]  +  out [ , "AX" ] 

#  - DOSEIN  =  Amount  absorbed  (mg) 

DOSEIN  <-  out [, "Alnh"] 

#  - MASSBAL  =  mass  balance  =  in  -  out  -  stored 

MASSBAL  <-  T MAS SO  +  DOSEIN  -  OUT  -  TMASS 
last  <-  MASSBAL [length (MASSBAL) ] 
names (last)  <-  "Mass  Balance  (mg)" 
print ( last ) 
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mb  <-  rbind (times,  rbind ( TMASSO ,  rbind (DOSEIN,  rbind(OUT,  rbind(TMASS, 
MASSBAL) ) ) ) ) 

#  View (mb) 


##################################################################### 
#  Plotting  routines 

source ( "plotrunWITHburnin . R" ) 
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APPENDIX  B.  SUPPORTING  R  SCRIPTS 


hypoxia2.R 

#  Modify  constants  to  simulate  hypoxia  scenario  2 

#  Changing  oxygen  in  atmosphere 

#  Start  with  fraction  02  (F02)  =  0.21. 

#  After  6  hour  burn  in,  reduce  F02  to  33%  or  0.07,  for  2  hours. 

#  At  8  hours,  return  to  baseline  F02  =  0.21. 

print ( "Sourcing  hypoxia2.R") 

#  Change  simulation  time,  oxygen  fractions,  and  change  times 

TSTOP  <-  12 

POINTS  <-  120 

CINT  <-  TSTOP  /  POINTS 

#  First,  turn  off  the  TCHNG  variable  (controlling  with  others) 
TCHNG  <-  Inf 

#  02  fractions 

Oxy  <-  c  ( 0 . 2 1 ,  0.07,  0.21,  0.21) 

#  Change  times  and  durations 

Dur  <-  c(6,  2,  6,  Inf) 

tau  <-  c(0,  Dur[l],  Dur[l]  +  Dur[2],  sum ( Dur [ 1 : 3 ] ) ) 


Fit.CB.R 


#  Fit.CB.R 

#  Get  relationship  between  CF  (free)  and  CB  (bound) . 
print ( "Sourcing  Fit.CB.R") 

Pa  <-  seq (0,150,5) 

source ( ' ~ / Documents/ 2017  Projects/02/ working/ Fit .Sa02 .orig.R' ) 
S  <-  function(p)  {  p**n/ (KD**n  +  p**n)  } 

caf  <-  function(p)  {  0.003*p  } 

cab  <-  function(p)  {  1 . 39*150*S (p)  } 

b  <-  cab ( Pa) 

e  <-  rnorm ( length (b) ,  0,  0.1) 
b  <-  b  +  e 
f  <-  caf (Pa) 

plot  ( f ,  b) 

bmaxO  <-  200 
kdO  <-  0.15 
nO  <-  2.5 

guess  <-  bmax0*f**n0/ (kd0**n0  +  f**n0) 
lines  ( f ,  guess ) 

#  Have  to  add  a  little  error  (e)  or  the  nls  algorithm  fails. 

#  So,  run  multiple  times  and  average  the  resulting  parameters, 
top  <-  100 

bsum  <-  0 
ksum  <-  0 
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nsum  <-  0 

for  (idx  in  l:top)  { 

fit  <-  nls (b  ~  Bmax*f**n/ (KD**n  +  f**n),  start=list (Bmax=bmaxO ,  KD=kdO, 
n=nO )  ) 

bsum  <-  bsum  +  coef(fit) [1] 
ksum  <-  ksum  +  coef(fit)  [2] 
nsum  <-  nsum  +  coef(fit)  [3] 

b  <-  cab (Pa) 

e  <-  rnorm ( length (b) ,  0,  0.1) 
b  <-  b  +  e 


summary (fit) 

pred  <-  predict (fit) 

lines  (  f ,  pred) 

Bmax  <-  bsum/top 
KD  <-  ksum/top 
n  <-  nsum/top 


plotrunWITHburnin.R 

#  calculate  burn-in  time 
end  <-  POINTS 

begin  <-  floor (POINTS*0 . 05) 

CT  <-  out [begin : end, "ATT" ] /VTT 
CTB  <-  out [begin : end, "ATB" ] /VTB 
CBr  <-  out [begin : end, "ABrT" ] /VBrT 
CBrB  <-  out [begin : end, "ABrB" ] /VBrB 
PctT  <-  out [begin : end, "ATT" ] /ATT0*100 
PctBr  <-  out [begin : end, "ABrT" ] /ABrT0*100 
t  <-  out [begin : end, "time" ] 
tau  <-  t  -  t  [  1 ] 

#plot(t,  PctT,  type="l",  xlab="Time  (hr)",  ylab="Blood  Percent  Oxygen") 
#plot(t,  PctBr,  type="l",  xlab="Time  (hr)",  ylab="Brain  Percent  Oxygen") 
plot  (tau,  out [begin : end, "A" ] ,  type="l",  xlab="Time  (hr)",  ylab="Arterial  Blood 
Concentration  (mg/L)") 

plot  (tau,  out [begin : end, "V" ] ,  type="l",  xlab="Time  (hr)",  ylab="Venous  Blood 
Concentration  (mg/L)") 

plot  (tau,  CT,  type="l",  xlab="Time  (hr)",  ylab="Body/Tissue  Oxygen 
Concentration  (mg/L)") 

plot  (tau,  CTB,  type="l",  xlab="Time  (hr)",  ylab=”Body/Tissue  Blood  Oxygen 
Concentration  (mg/L)") 

plot  (tau,  CBr,  type="l",  xlab="Time  (hr)",  ylab=”Brain  Oxygen  Concentration 
(mg/L) " ) 

plot (tau,  CBrB,  type="l",  xlab="Time  (hr)",  ylab="Brain  Blood  Oxygen 
Concentration  (mg/L)") 

plot  (tau,  out [begin : end, "HIF" ] ,  type="l",  ylab="HIF") 
plot  (tau,  out [begin : end, "HIFN" ] ,  type="l",  ylab="HIFN") 
plot (tau,  out [begin : end, "TC" ] ,  type="l",  ylab="TC") 
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APPENDIX  C.  MODEL  PARAMETERS 


Qc 

Qp 

QBr 

Ca 

Cf 

Cb 

Cv 

Cf 

Ci 

Cx 

CvBr 

Vc 

Vn 

A 


ABr,T,  CBr.T 

VMax,  KM 

VMaXBr,  Vresp,  Kresp 

PA 

PBr 

PB 

Bmax,  KD,  n 

Hgb 

Pa02 

Sa02 

Kphd 

KmHIF,  Km02 


Khif 
T max.  Kt 


Kndeg 

Kbm 

Ktc 

F 


Cardiac  output  (total  blood  flow) 

Pulmonary  blood  flow 
Brain  blood  flow 

Arterial  blood  oxygen  concentration  (=Cf  +  Cb) 

Free  oxygen  concentration  in  blood 
Bound  oxygen  concentration  in  blood 
Venous  blood  O2  concentration 
Free  oxygen  concentration  in  blood 
Inhaled  air  O2  concentration  (alveolar) 

Exhaled  air  O2  concentration 
O2  blood  (venous)  concentration  leaving  brain 
Brain  intracellular  compartment  volume 
Brain  nuclear  compartment  volume 
Amount  of  O2  in  arterial  blood 

Amount  and  concentration  of  O2  in  brain  intracellular  compartment 

Michaelis-Menten  parameters  for  O2  metabolism 

Michaelis-Menten  parameters  for  O2  metabolism  in  brain 

Permeability- surface  area  product  for  O2  uptake  into  brain 

Blood-brain  partition  coefficient  for  O2 

Air/blood  oxygen  partition  coefficient 

Hill  coefficients  for  O2  hemoglobin  binding 

Amount  of  hemoglobin 

Partial  pressure  of  O2  in  blood 

Fraction  saturation  of  hemoglobin 

PHD-dependent  degradation  rate  constant  for  HIF-la  in  the  cytosol 

HIF-la  and  O2  half- saturation  constants  for  PHD-dependent  degradation 

rate  constant  for  HIF-la  in  the  cytosol 

Zero  order  synthesis  rate  of  HIF-la  in  the  cytosol 

Maximum  rate  and  half-saturation  concentration  for  transport  of  HIF-la 

from  cytosol  to  nucleus 

Degradation  rate  of  HIF-la  in  the  nucleus 

Rate  constant  for  binding  of  HIF-la  with  hypoxia  responsive  element 
(HRE)  on  the  DNA 

Spontaneous  decay  rate  of  the  transcription  signal  TC 
Oxygen  unit  conversion  factor  (nM/s  to  mg/hr) 
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LIST  OF  ACRONYMS 


BSA 

bovine  serum  albumin 

CB 

carotid  body 

CBP 

cAMP  response  element  binding  (CREB)  protein 

COHb 

carboxyhaemoglobin 

2,3-DPG 

2,3-diphosphoglycerate 

ELISA 

enzyme-linked  immunosorbent  assay 

FIH 

factor  inhibiting  hypoxia  inducing  factor 

Hb 

hemaglobin 

HIF 

hypoxia  inducing  factor 

HIFN 

hypoxia  inducing  factor  in  the  nucleus 

HRE 

hypoxia  responsive  element 

IHC 

Immunohistochemical  (staining) 

OBOGS 

on-board  oxygen  generation  system 

ODE 

ordinary  differential  equation 

PA 

permeability  area  cross  product 

PBPK 

physiologically-based  pharmacokinetic  (model) 

PHD 

prolyl  hydroxylase 

p02 

partial  pressure  of  oxygen  (mmHg) 

ss 

steady-state 

TAD 

transactivation  domain 

TC 

transcription  signal 

VHL 

von  Hippel-Lindau  tumor  suppressor  protein 

VEGF 

vascular  endothelial  growth  factor 

WMH 

white  matter  hyperintensities 
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